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1.  Introduction 

One  of  the  most  promising  approach  for  remote  sensing  is  polarimetry.  Polarimetry 
assumes  the  exposure  of  the  medium  with  electromagnetic  radiation  with  a  given  polarization 
and  the  subsequent  study  of  the  transformation  of  this  polarization  as  a  result  of  interaction 
between  incident  radiation  and  studied  medium. 

It  is  known  that  if  the  interaction  of  radiation  with  the  medium  is  linear,  then  it 
mathematically  can  be  represented  by  a  system  of  linear  equations: 

Eout  =  TEin  (1) 

where  E,n(out)  -  vectors  representing  the  polarization  of  input  and  output  radiation;  T  -  is  a 
matrix  that  defines  a  linear  transformation  of  polarization  state  of  input  radiation  into  the  output 
one. 

When  one  intends  to  find  the  unknown  parameters  of  the  output  radiation  E'n<out)  for  a 
given  T ,  this  is  the  so-called  “direct  problem”.  We  interested  in  the  “inverse  problem”,  the 
essence  of  which  is  finding  the  explicit  form  of  the  matrix  T .  It  is  clear  that  the  matrix  T 
contains  all  information  about  the  properties  of  the  medium  in  implicit  form.  To  extract  this 
information  a  number  of  relevant  matrix  models  have  been  developed.  The  most  perspective  are 
the  models  based  on  the  singular  and  the  polar  decompositions  [1-5].  Using  the  singular  and  the 
polar  decomposition  with  respect  to  the  operator  T  gives  the  possibility  to  represent  an 
investigated  media  as  a  finite  set  of  layers  with  well-studied  and  easily  interpretable  anisotropic 
properties.  In  the  case  of  singular  value  decomposition  the  medium  is  represented  as  a  set  of  four 
layers  -  one  with  a  linear  dichroism,  one  with  the  optical  activity,  and  two  layers  with  linear 
birefringence  [6].  By  polar  decomposition  the  same  medium  is  represented  as  two  layers,  which 
are  the  elliptical  polarizer  and  an  elliptical  phase  plate  of  general  forms  [2-5],  Thus,  the  complex 
anisotropy  of  the  arbitrary  deterministic  medium  is  decomposed  into  simple  components 
describing  separately  the  transformation  of  phase  and  amplitude  of  the  input  radiation  by  studied 
medium. 

The  fact  that  the  model  basing  on  singular  value  decomposition  contains  two  components 
with  the  same  types  of  anisotropy  (linear  birefringence)  complicates  the  unique  interpretation  of 
decomposition  results  and  thus  leads  to  a  relatively  rare  using  of  this  model.  Polar  decomposition 
is  used  much  more  frequently,  see  for  example  Refs.  4,5,7  and  9.  However,  from  our  point  of 
view,  the  parameters  which  are  obtained  in  polar  decomposition,  i.e.,  diatenuation  vector, 
polarizance  and  retardation,  have  also  no  clear  physical  interpretation. 

In  scope  of  this  project  we  used  a  model  based  on  the  so-called  generalized  equivalence 
theorem  [10].  In  concordance  with  this  theorem  arbitrary  deterministic  medium  can  be  presented 
as  a  product  of  four  basic  types  of  anisotropy:  linear  phase  and  amplitude  and  circular  phase  and 
amplitude.  Physics  of  these  types  of  anisotropy  is  studied  well  and  parameters  characterizing 
these  type  of  anisotropy  have  clear. 

Note  that  all  mentioned  multiplicative  layered  models  of  anisotropic  media  consist  of  the 
components  which  characterized  by  orthogonal  eigenpolarizations.  However  the 
eigenpolarizations  of  the  initial  medium  are  not  necessarily  orthogonal.  Elliptical  polarizer  in 
polar  decomposition  also  has  orthogonal  eigenpolarizations.  At  the  same  time,  the  medium 
represented  by  a  sequence  of  layers  with  dichotic  properties:  linear  and  circular  dichroism,  in 
general  has  a  non-orthogonal  eigenpolarizations.  The  absence  of  dichroism  in  medium  results 
immediately  in  orthogonality  of  eigenpolarizations. 

The  cases  of  orthogonal  eigenpolarizations  usually  receives  increased  attention  in  the 
literature,  and  it  has  been  a  milestone  for  classifying  polarization  elements  and  studying  their 
properties  [1,2,6,11],  However,  in  scope  of  multiplicative  matrix  models  there  is  no  systematic 
study  of  the  conditions  under  which  the  arbitrary  anisotropic  medium  will  have  orthogonal 
eigenpolarizations  in  general  case.  Derivation  of  these  conditions  will  determine  the  structure 
and  symmetry  of  the  matrix  model  of  corresponding  class  of  media. 
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Thus,  the  main  goal  of  this  study  is  derivation  of  the  conditions,  under  which  arbitrary 
crystalline  anisotropic  medium  has  in  general  the  orthogonal  eigenpolarizations.  In  addition  we 
(i)  developed  the  scheme  of  polarimeter,  which  is  optimally  fit  for  the  measurement  of  structures 
of  Mueller  matrices  with  orthogonal  eigenpolarizations;  and  (ii)  assembled  the  breadboard  of  ID 
polarimeter,  which  meets  the  requirements  derived  in  scope  of  this  research. 

Because  of  unknown  structure  of  the  matrix  model  for  media  with  orthogonal 
eigenpolarizations  in  general  case  developed  polarimeter  should  provide  an  “operation 
flexibility”  between  various  measurement  strategies  for  optimal  measurement  of  given  structures 
of  matrices.  One  way  to  do  this  is  using  in  input  and  output  channels  of  polarimeter  the 
polarization  transducers  with  universally  controllable  parameters. 

Currently  a  number  of  schemes  to  measure  the  Stokes  parameters  of  electromagnetic 
radiation  (Stokes-polarimeter)  have  been  proposed  and  implemented  [11-15],  Key  elements  for 
these  systems  are  polarization  transducers  with  controllable  parameters.  Basing  on  approach  to 
control  of  polarization  parameters  the  Stokes-polarimeters  can  be  divided  into  two  classes: 
mechanically  controllable  and  electrically  controllable.  Mechanically  controlled  transducers  are 
the  ones  which  parameters  are  mechanically  altered  (moved  or  rotated  polarizers,  phase  plates, 
etc.).  Schemes  of  Stokes-polarimeter  with  mechanically  controlled  polarization  transducers  are 
widely-spread  due  to  simplicity  in  implementation  and  adjusting.  However,  they  are 
characterized  by  significant  disadvantages  due  to  impossibility  of  providing  the  necessary  level 
of  accuracy  and  performance. 

Among  the  electrically  controlled  transducers  it  can  be  pointed  out  electro-,  magneto-, 
acousto-optical  cells  etc.  Anisotropic  properties  of  these  transducers  are  changed  by  applying  an 
external  electric  (magnetic)  field  which  can  be  accomplished  with  a  quite  high  speed  and 
accuracy  without  mechanical  moving.  Disadvantages  of  electrically  controlled  polarization 
transducers  are:  1)  nonlinear  dependence  of  parametres  on  the  applied  external  field,  and  2) 
limitation  of  the  range  of  anisotropy  values,  3)  as  a  rule  only  one  parameter  can  be  controlled 
(the  orientation  axis  of  anisotropy,  values  of  anisotropy  etc.).  These  disadvantages  somewhat 
complicate  the  adjustment  of  Stokes-polarimeter  and  require  the  using  of  more  than  one 
transducer,  or  changing  their  orientations  to  expand  the  range  of  polarization  transformation. 
Recently  transducers  on  liquid  crystals  (LC)  cells  [16-29]  become  popular.  LC  cell  consists  of  a 
layer  of  liquid  crystal  located  between  two  transparent  electrodes  controlled  by  external  voltage. 
This  cell  exhibits  birefringence  which  axis  is  normal  to  the  propagation  direction.  Thus,  the  LC 
cell  is  a  wave  plate  (LC  retarder)  which  phase  shift  (birefingence)  is  controlled  by  voltage. 
Coupling  the  relatively  low  price  and  satisfactory  performance  the  LC  cell  allows  to  get  a  wide 
range  of  changes  of  birefringence  by  controlling  voltage  near  15-30V.  Two  consequently  placed 
LC  cells  provide  a  range  of  polarization  modulation  of  radiation  enables  to  determine  all  four 
Stokes  parameters. 

2.  Methods,  Assumptions  and  Procedures 

2.1  Jones  and  Mueller  matrix  methods.  Spectral  problem 

To  describe  the  linear  interaction  of  polarized  radiation  with  the  medium,  see  Eq.(l),  the 
Jones  and  Mueller  matrix  methods,  which  is  uniquely  related  in  case  of  a  homogeneous 
anisotropic  media  [30],  are  used. 

When  Jones  matrix  method  is  used  than  in  Eq.(l)  E,n<out)  is  the  Jones  vector  of  input 
(output)  radiation,  and  T  denotes  the  Jones  matrix. 

Jones  matrix  T  (2x2  matrix  with  complex  elements  tmn)  describes  anisotropic  properties  of 
homogeneous  medium: 


T  = 


21 


Where  Kn  =\tmn  leXPH^J' 


(2) 
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Eigenpolarizations  x  °f  such  matrix  can  be  obtained  as: 


X\,2 


1  ?22  hi  —  "\/(^22  ^  +  4tj2  ’  ?2i 

2  T1X  ; 


(3) 


where  /  =  Ex  /  Ey  -  complex  variable  [30];  E  -  components  of  the  Jones  vector  E  . 
For  orthogonal  eigenpolarizations  the  following  relation: 


Z1Z2  =  -1 


(4) 


is  satisfied. 

In  experimental  studies  the  Mueller  method  is  used  for  description  of  interaction  between 
electromagnetic  radiation  and  medium  because  it  operates  with  intensities  of  radiation  that  can 
be  directly  measured. 

When,  in  scope  of  the  Muller  matrix  method  Eq.(l)  can  be  rewritten  as: 

Sout  =  M  Sin  (5) 


where  Sin(out)  -  denotes  input  (output)  Stokes  vector. 

The  definition  of  Stokes  vector  is  follows  : 

/ 

Ip  cos(2/?)  cos(2£-) 
Jpsin(2/?)  cos(2£-) 
Ip  sin(2£-) 


(6) 


where  /  -  overall  intensity  of  radiation;  p  -  polarization  degree;  3  -  azimuth  and  e  -  ellipticity 
angle  of  polarization  ellipse. 

It  can  be  seen  from  Eq.(6)  that  the  elements  of  Stokes  vector  (called  the  Stokes  parameters) 
have  dimensions  of  intensity.  Thus,  the  Stokes  parameters  and  Mueller  matrix  elements  can  be 
measured  directly  in  experiment.  Besides,  the  Stokes  parameters  describe  either  completely 
polarized  (p=  1),  partially  polarized  (0<p<l)  and  depolarized  (p=0)  radiation.  Similarly,  the 
Mueller  matrix  can  represent  both  homogeneous  and  inhomogeneous  media. 

In  accordance  with  Eq.(5)  the  Mueller  matrix  M  is  a  4x4  matrix  with  real  elements 


m  11 

mn 

mu 

ml4 

m2 1 

m2  2 

m23 

m24 

m3l 

m32 

m33 

m34 

mAX 

m43 

m44 

(7) 


This  matrix,  as  well  as  Jones  matrix,  describes  completely  anisotropic  properties  of 
homogeneous  medium  for  a  given  input  and  output  (scattering)  directions  and  wavelength  of 
input  radiation. 

Direct  solving  of  the  spectral  problem  [31]  in  scope  of  Mueller  formalism,  i.e.,  finding  the 
conditions  on  Mueller  matrix  elements  for  eigenpolarizations  to  be  orthogonal,  is  quite 
complicated  task  because  of  Mueller  matrix  dimension.  However,  this  problem,  as  it  is 
demonstrated  below,  can  be  solved  in  scope  of  the  Jones  formalism  both  in  terms  of  matrix 
elements  and  in  terms  of  anisotropy  parameters  (introduced  below,  see  subsection  2.3).  Taking 
into  account  the  fact  that  Mueller  M  and  Jones  T  matrices  for  homogeneous  medium  are 
interconnected  by  relation 
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(8) 


M  =  A(T  0  T*  )A_1 , 

where  *  -  conjugation;  0  -  Kronecker  product;  and 


10  0  1 

10  0-1 
0  110 
0  i  -  i  0 


all  results,  which  are  obtained  for  the  Jones  formalism,  can  be  translated  to  the  Mueller 
formalism  with  Eq.(8).  Note  that  the  main  condition  for  that  is  the  medium  under  consideration 
does  not  depolarize  input  radiation. 

It  is  important  to  note  that  any  2x2  matrix  could  be  called  “Jones  matrix”,  i.e.,  arbitrary 
2x2  matrix  with  complex  elements  describes  always  a  physical  realizable  transformation  of 
polarization  trough  the  Eq.(l).  The  same  can  not  to  be  said  about  any  4x4  matrix  with  real 
elements.  To  be  named  as  “Mueller  matrix”  this  matrix  has  to  meet  an  ample  of  requirements 
[32]. 

Due  to  above  conditions  the  solving  of  spectral  problem  even  for  Mueller- Jones  matrix 
(which  is  determined  by  Eq.(8))  case  is  generally  difficult.  We  can  write  it  in  the  following  form: 


1 

0 

0 

0 


(9) 


2.2  Generalized  equivalence  theorem 

For  analysis  of  anisotropic  parameters  of  medium  we  use  the  so-called  generalized  matrix 
equivalence  theorem  [33].  According  to  this  theorem  an  arbitrary  Jones  (Mueller- Jones)  matrix 
of  crystalline  anisotropic  medium  can  be  represented  as  the  product  of  the  Jones  matrices  of  four 
basic  types  of  anisotropy:  linear  dichroism  (linear  amplitude  anisotropy)  TLA ,  circular  dichroism 
(circular  amplitude  anisotropy)  TCA ,  linear  birefringence  (linear  phase  anisotropy)  TLP ,  and 
optical  activity  or  circular  birefringence  (circular  phase  anisotropy)  TCP : 


T  = 


rp  CP  rp  IP  rp  CA  rp  LA 


(10) 


The  Jones  and  Mueller  matrices  of  basic  types  of  anisotropy  are  well-known  [30]  and 
presented  below. 

Jones  matrices: 


cos 2  (y)  +  sin 2  (y)P  cos(y)  sin(y)(l  -  P ) 
cos(y)  sin(y)(l  -  P)  sin 2  (y)  +  cos 2  (y  )P 
1  -i-R 
i-R  1  ’ 


T 


LP 


T 


CP 


cos2  (a)  +  sin2  (a)e”'A  cos(a)  sin(a)(l  -  e”!A ) 
cos(a)  sin(a)(l  -  e”'A )  sin2  (a)  +  cos2  (a)e_,A 
cos((p)  sin(cp) 

-  sin((p)  cos((p)_ 


(11) 


Mueller  matrices: 
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1 +  R2  0  0  ±2  R 

0  1  -R2  0  0 

0  0  1  -R2  0 

±2  R  0  0  1 +  R2 

10  0  0 
0  cos(2  <p)  sm(2<p)  0 

0  -  sin(2^>)  cos(2#>)  0 

0  0  0  1 

(12) 

1 +  P2  (l-.P2)cos(2y)  (1-P2)sin(2y)  0 

(\-P2) cosily)  cos2(ly\\  + P2)  +  lsm2(ly)P  cos(2y)sm(ly\\- P)2  0 

(1-P2)sin(2y)  cos(2y)sin(2y)(l-.P)2  sin  2(2y)(l  +  P2)  +  2  cos2  (2 y)P  0  ’ 

0  0  0  IP 

10  0  0 
0  cos2(2a)+  sin2(2a)cos(A)  cos(2a)sin(2a)(l  -  cos(A))  -  sin(2a:)sin(A) 

0  cos(2<z)sin(2<z)(l  -  cos(A))  sin2(2<z)  +  cos2(2<z)cos(A)  cos(2«)sin(A) 

0  sin(2<z)sin(A)  -cos(2«)sin(A)  cos(A) 

where:  value  of  linear  dichroism  is  in  range  P40;l]  and  the  azimuth  of  maximum  transition  is  in 
range  7^-tt/2;tt/2];  value  of  circular  dichroism  is  in  range  R  <|-1;1];  value  of  linear  birefringence 
is  in  range  A  ^0;2tt]  with  it’s  fast  axis  orientation  in  range  a^-7t/2;7t/2];  value  of  optical  activity 
is  in  range  ^-7t;7i]. 

Parameters,  which  characterize  basic  types  of  anisotropy,  have  clear  physical  meanings. 
They  associate  with  time  and  spatial  non-locality  of  the  medium  response  on  input  light,  i.e., 
with  time  and  spatial  dispersion.  Last  ones  determine  the  form  of  Maxwell’s  constitutive 
relations.  Thus,  generalized  polarimetric  equivalence  theorem  combines  both  mathematical 
generality  and  physical  interpretability. 

2.3  Symmetry  of  matrix  model  of  medium  with  orthogonal  eigenpolarizations. 

2.3.1  Conditions  on  matrices  elements 

Eigenpolarizations  are  those  polarization  states  of  light  that  do  not  change  when  passing 
through  a  medium.  The  amplitude  and  the  overall  phase  of  the  beam  of  light  with  an 
eigenpolarization  do,  however,  change.  These  changes  are  described  by  the  corresponding 
eigenvalues.  In  optics  and  electrodynamics  the  crystalline  medium  is  characterized  by  the  types 
of  eigenpolarizations  that  this  medium  possesses.  Because  of  that,  ascertainment  of  generalized 
conditions  for  orthogonality  of  medium’s  eigenpolarizations  allows  determining  the  structure 
and  symmetry  of  the  matrix  model  for  such  class  of  media. 

Using  Eq.(3)  we  have 

X\Xi  =  ~~~  ■  (13) 

1 21 

Combining  orthogonality  condition  in  terms  of  the  complex  variables  Eq.(4)  with  Eq.(3) 
we  get  first  relation  between  non-diagonal  elements  of  the  Jones  matrix  in  general  case: 


8 
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(14) 


tX2  t2l 


Next,  it  can  be  seen  that: 


Izi + z2f  =  |;&  f + 1*2  f  +  *1*2* + *2*i*  =  I*  f  +r\~2  = 

M 

i  Y 

w-rh 

m) 


k.r-02 


*i 


=  m\ 


If  we  rewrite  the  relation  for  eigenvectors  as: 


*1,2  =  » where  e,  =t22 -tu;  e2  =^22  “hi)2  +4M2i  > 

zr2i 


then  Eq.(15)  takes  the  form: 


el 

e,  +e2 

<N 

1 

^21 

2t2l 

2^21 

or 


h  +ei  I 

\ei  e2 

K 

=  2. 


Let  us  present  the  parameters  e\  and  e 2  in  the  exponential  form: 

ei,2  =  A, 2  exp(-/>12  ). 

Then,  after  some  mathematics,  from  Eq.(17.b)  we  obtain: 

cos(^!  -y/2)2  =1, 

¥1  =  ¥2  ±nn . 

Combining  the  equations  Eq.(16)  and  Eq.(19.b),  for  the  case  n  =  0  we  can  write: 
\e2 1 2  exp(2/  Vi ) = \l2i  -  hi  |2  exp(2/ ¥1 )+ 4|h2  |h2 1  expK^i2  +  K )]  • 


(15) 


(16) 


(17.a) 

(17.b) 

(18) 

(19. a) 
(19  b) 


(20) 


From  Eq.(20)  it  can  be  deduced  that  2  ^1=^  12+^21  ±  nn  ox,  taking  into  account  Eq.(16),  it 
is  transformed  into 


2022-ll  ^12  ^21  — 


(21) 


For  sum  of  matrix  elements  C 1  and  tu,  allowing  for  Eq.(21),  we  can  write: 
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,|exp(/^12+21)  =  |/12|[(cos^12  +cos^21)+/(sin^ 


^21  +  ^12  |^21  +  hi 

:  2|f12 1  cos  expf  i  ' 

2  y  2 


>12  +  Sln  (p2\  . 


(22) 


y 


from  which  it  follows: 


24+12  4  4  — 


(23) 


excluding  the  case  when  *21  =  - 1\2- 

Taking  into  account  Eq.(14),  one  more  relation  on  phases  of  difference  of  matrix  elements 
*2i  and  *12  can  be  also  written: 


i(*2i  -  *12 )  =  |/(/2i  - 1\2 1  exp0>2l_l2 )  =  i\tn  |[(cos  </>n  -  cos  </>2l )  +  /(sin  (j)n  -  sin  (j>2[ )]  = 


=  2 \tn  |  sin 


4  4 


exp 


•  4  4 


(24) 


From  Eq.(24)  it  follows: 


2$(21-12)  =4  +  4 


(25) 


Using  interrelation  Eq.(8)  between  Jones  and  Mueller  methods  and  description  of  spectral 
problem  as  Eq.(9)  we  can  also  study  symmetry  of  the  Mueller  matrix  of  medium  with  orthogonal 
polarizations. 

In  particular,  from  Eq.(8)  it  follows: 


r  21  — 


mn  —m2 1  +ml2  -m22 


’  1^12  — 


mu  +  m2l  -  ml2  -  m22 


hi  — 


mu  +  m2l  +  ml2  +  m22 


(26) 


’  r  22 


mu  -  m2l  -  mn  +  m22 


cos(</>l2)  = 


mU  +  m23 


2 hi  r 


12 


cos(^21)  = 


m31  +  m32 


sin(4)  = 


ml4  +  m24 


2 hi  r 


12 


2knk 


sin(^21)  = 


m4l  +  m42 


cos(^22)  = 


21 1 

m33  +  m44 


2kllk: 


(27) 


2knk 


sin(^22)  = 


ill 


m43  -  m34 


22 


2knk 


12 


Here  we  assumed  that  all  phases  of  Jones  matrix  elements  Eq.2  are  normalized  on  phase  of 
the  first  matrix  element  (i.e.  (j)mn  — »  (/)mn  ~</>n)-  Thus,  in  this  case  the  phase  of  element  mn  is 
->  0- 

Condition  Eq.(14)  for  orthogonality  of  eigenpolarizations  of  medium  in  terms  of  Jones 
matrix  elements  can  be  transformed  using  Eq.(24)  into  relation: 

mu  -  m2l  +  ml2  -  m22  =  mn  +  m2l  -  ml2  -  m22 , 


or 
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mu  =  m2l 


(28) 


Condition  Eq.(21)  in  terms  of  Jones  matrix  elements  can  be  transformed  using  Eq.(27)  into 
relations: 


*mn  =  Kn  \{™s((/)mn )  +  i  sin(^m„ )) ; 


tg(^22-l  1  )  — 


sin(4  +^21)  _ 
COS(^12  +  (f>2\ ) 

Im(t22  —  tl  j ) 


tg 


2  arctg 


Im(t22  ( i ) 


Re(t22  t\  i ) 


m43  -  m34 


72 


tn)  m33  +  m44  -  mn  -  m2l  -  ml2  -  m22 


R  e(t22 

cos(^i2  +  </>2l )  =  cos(^i2 )  cos(4i )  -  sin(^i2 )  sin(^21 ) ; 

sin(^i2  +  (j)2X )  =  cos(^i2 )  sin(^2i )  +  cos(^2i )  sin(^i2 ) ; 

C0S(^12  +  ^  )  =  m!3  +  m23  ffl31  +ffl32  +  ™14  +  W24  ™41  +  «42 


2|h 1 1|^12 1  ^11 1|^21 1  2|tn|ti2|  2|tu|t2i| 


sin(^2  +^21)  = 


W13  +  "*23  "Cl  +  "*42 


W14  +  m2A  "Cl  +  "*32 


||^i2 1  2|Ci||Ci|  2|tn|ti2|  2|tn|t2i| 

or  in  terms  of  Mueller  matrix  elements 

("Ci  +  m42)(m23  +mn)~  (ml4  +  m24)(m32  +  m3l ) 


("C 1  +  m42  )  ("*14  +  m24  )  +  ("*23  +  "*13  )  ("*32  +  "*31  ) 


=  tg 


2  arctg 


"C3  -m34 


^  mu  +  m2 1  +  m ,  2  +  m22 


n  (29) 


m33  ~  m 


44  J 


Expression  Eq.(29)  is  evidently  inconvenient  for  practical  use.  Thus,  we  propose  to  use 
relations  2$1+22  =2^2|+|2)  =  2$(21+l2)  =<fr2  +  ^2I  ±nn  (see  Eq.(21),  Eq.(23),  Eq.(25)).  Then: 

Im(;(t21  ~  tl2 ))  ^  Im(t2i  +  tj2 ) 

Rc(/(/2|  — 1[2))  Re(t2i  +h2) 


*(Cl  t j2  )  —  tj2 


m3\  +  m32  "*13  +  "*23 


m4 1  +  m42  +  m24 


2|hi|h2|  2|tn|ti2|  2|tn|t12|  2|?u  ||?i2 1 


=  -r— 3 1  +i 

2|hi| 


^21  +  ^12  kl2 


mn  +  m23  |  m31  +  m32  | 


?/  h  ?k  k 

^rn  ri2 


12 


^14  +  m24  ^41  +  m42 


2  k  Jk 


ii  r  12 


2kn  k 


=  ^  [m31  +  m32  +  m13  +  m23  +  /(m41  + 

ii  | 

Im(/(^21  -tl2))=  Im(f21  +  f12) 
Re(/(^2i  “^12))  R^(^2i  ^12 ) 


^42  “  ^14  “  ^24 


12 1  y 

]; 
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or 


mX3+m23  -m3l  -m32 


m4\  +  m42  ~  mi4  ~  m24 


m4i  +  m42  +  mu  +  m24 


mn  +  m23  +  m3l  +  m32 


(mi3  +  m22  f  ~  (mn  +  m32  f  =  (m4\  +  m42  f  “  (™14  “  W24  )" 


(30) 

(30') 


Finally  we  have 


W\  3  +  ^23 


-  m 


31 


m 


32 


m4i  +  m42 


■nt 


14 


■m 


24 


m 


43 


m 


34 


m4i  +  m42  +  m14  +  m24 


/77 , 3  +  m23  +  /773 ,  +  m32 


mn  +  m21  +  m\2  +  m22  ~  m33  ~  m44 


(31) 


Spectral  problem  in  scope  of  Mueller  formalism  can  be  expressed  in  form  Eq.(9)  where  the 
fact  that  eigenpolarizations  S  of  the  Mueller  matrix  are  orthogonal  was  already  taken  into 
account.  In  other  words  Eq.(9)  holds  for  following  normalized  Stokes  vectors  : 


S,= 


rn 

f  1  3 

—s~ 

2 

II 

(N 

Xfl 

2 

^3 

~S3 

v  S  A  j 

<  j 

(32) 


For  Stokes  vectors  Eq.(32)  one  can  write: 


*12 


m3\ 

\m4X 


m 


32 


m 


*'14 


m2 ,  m22  m23  m24 


33 


m 


34 


m42  m43  m 


s. 


44  J  \S4  J 


(  mn  +  mus2  +  mus3  +  ml4s4  ^ 
m2l  +  m22s2  +  m23s3  +  m24s4 
in  3 ,  +  m32s2  +  m33s3  +  m34s4 
vm41  +  m42s2  +  m43s3  +  m44s4  j 


=  (mu  +  m12s2  +  m13s3  +  m14s4) 


^2 

^3 


m . 


/77 


/77 


31 


Vm41 


777 


12 


21 


772 


777 


777 


777 


13 


'22 


777 


32 


777 


33 


'42 


777 


'43 


(33) 


m 


14 


23 


777 


'24 


777 


34 


777 


1 

“^2 


A  f  ™n-ml2s2-ml3s3-ml4s4^ 


44  y  “’4  J 


m 


21 


-m23s3-m24s4 


m3\  ~  m32S2  ~  m32S3  ~  m 34s 4 


V«41  -  m42S2  ~  V3  ~  m44S4  j 


=  (mu-mns2  -ml3s3-ml4s4) 


1 

~S2 

-s. 


~S4J 


From  Eq.(33)  it  follows: 


^  mu+  ml2s2  +  ml3s3  +  ml4s4  ' 

m2l  +  m22  S2  +  m23S3  +  m24S4 

m3 1  +  m32.v2  +  m33s3  +  m34s4 
Km4l  +  m42s2  +  m43s3  +  m44s4  , 


(mn  +ot1252  +  mus3  +ml4s4) 


rn 

fOl 

*2 

0 

S3 

0 

^  $A  j 

,0, 

(34) 


r  mn  -  mus2  -  ml3s3  -  mus4  ^ 


W31  ~m32S2  ~m33S3  ~m34S4 


V«4i  -  m42s2  -m43s3-m44s4j 


~{mn  -ml2s2  -ml3s3  -  mus4 ) 


f  1  > 

^0^ 

~S2 

0 

S3 

0 

,0, 

(35) 


or 
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m2i 

+ 

m22S2 

+ 

m23s3 

+ 

m24s4  - 

-s2(mu 

+ 

m\2S2 

+  mn 

s3 

+  mu 

sa) 

=  0 

m3l 

+ 

m32S2 

+ 

m33s3 

+ 

m34s4  - 

~s3  On 

+ 

m\2S2 

+  ml3 

s3 

+  mu 

=  0 

m4l 

+ 

m42s2 

+ 

m43s3 

+ 

m44s4  - 

-s4(mu 

+ 

m\2S2 

+  ml3 

s3 

+  ml4 

sa) 

=  0 

m2l 

- 

m22S2 

- 

m23s3 

- 

m24s4  - 

-s2(mu 

- 

m\2S2 

M\  3 

s3 

~mU 

*4) 

=  0 

m3 1 

- 

m32S2 

- 

m33s3 

- 

m34s4  - 

-s3(mu 

- 

m\2S2 

—  ml3, 

s3 

-mw 

S4) 

=  0 

m4l 

- 

m42S2 

- 

m43s3 

- 

m44s4  - 

~s4(mu 

- 

m\2S2 

~mu 

s3 

-mu 

sa) 

=  0 

By  adding  Eq.(34)  to  Eq.(35)  we  obtain: 

m2l  =  s2(mns2  +  mus3  +  ml4s4), 

m3l  =  s3  ( mns2  +  ml3s3  +  mu  s4 ),  (36) 

m4l  =s4(ml2s2  +  ml3s3  +ml4s4). 

Multiplying  Eq.(36)  by  S2,  s3,  S4,  respectively,  and,  taking  into  account  that  for  completely 

(r\  r\  \  1/2 

s2  +  s3  +  s4  1  =  1 ,  we 

have: 

s2  (m2l  ~ml2)  +  s3  (m3l  -  ml3 )  +  .s4  (m41  -  ml4  )  =  0 ,  (37) 


or 


m2l  {m2l  ~mn)  +  m2l{ 


mn  +  mn  +  m\4 

mi2m2l  +  mi3m31  +  m\4m4l 


(m4l-ml4)  =  0, 

(38) 

=  1. 

(39) 

Taking  into  account  that  for  Mueller- Jones  matrix  the  following  condition 


ml22  +  ml32  +  ml42  =  m22  +  m32  +  m42 , 


(40) 


is  true  [34],  from  Eq.(39)  and  Eq.(40)  we  can  write  that: 


2  2  2  2  2  2 

mn  +  mn  +  m\4  +  m2\  +  m2\  +  m4\  =  2 (ml2m2l  +  mX3m3X  +  ml4m4l )  : 

\2 


(mn  -  m2l )  +(ml3  -m3l)  +  (ml4  -  m4lf  =  0 


(41) 


from  Eq.  (4 1 )  it  results  that : 

mu  —  tw2i  , 
/77 J 3  =  m3l , 

mu  =m4V 


(42) 


From  Eq.(42)  it  follows  that  for  complete  description  of  anisotropy  of  medium  with 
orthogonal  eigenpolarizations  the  knowledge  of  second,  third  and  fourth  columns(rows)  of 
Mueller  matrix  is  sufficient.  This  result  can  be  used  during  optimization  of  polarimeter  for 
studying  of  given  polarization  class  of  media. 
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2.3.2  Orthogonality  conditions  in  terms  of  anisotropy  parameters. 

To  derive  orthogonality  condition  in  terms  of  anisotropy  parameters  we  use  the  mentioned 
matrix  model  of  arbitrary  homogeneous  anisotropy  (deterministic)  Eq.(10),  that  has  recently 
been  presented  in  [33]. 

The  relations  Eqs.(16),  Eq.(22)  and  Eq.(24)  appear  useful  to  rewrite  in  the  form: 


t22  -tn  =  (ax  +ibl)ex p 


f  .ts> 

-l  — 

.  2y 


(43) 


t2 1  -  ti2  =  ( a2  +  ib2  )exp 


-i  — 

v  2  j 


(44) 


t2l  +tl2  =  (a3  +  ib3)exp 


f  aA 


(45) 


V  ± J 


Where,  in  compliance  with  Eqs.(10)  and  Eq.(ll),  parameters  ai  and  bt  have  the  following 


form: 


ax  =  (l  +  T,)f?sinysin(2a-(p)-(l-T’)cosYCOs(2y-cp), 
bx  =  (l  -  P)R  cos  y  sin(2y  -  (p)  -  (l  +  P)sin  cos(2a  —  q>), 
a2  =  (l-T,)f?sinyCOs(2a-2y-cp)-(l  +  T,)cosYsin(cp), 


b2  =  (l  +  />)f?cosyCOs(cp)-(l-T,)sinysin(2a-2y-cp), 
a3  =  (l  +  P)R  sin  y  cos(2a  -  tp)  +  (l  -  P)cos  y  sin(2y  -  (p) , 
b3  =  (l  -  P)R  cos  y  cos(2y  -  (p)  +  (l  +  />)sin  y  sin  (2a  -  cp) . 


Combining  Eqs.(44)  and  (45),  equation  (14)  can  be  written: 


(46) 


a2a3  +b2b3  =  0. 


(47) 


As  it  follows  from  Eq.(21),  Eq.(23)  and  Eq.(25) 

^22-11  =  ^12+21  =  ^/(21-12)  > 


(48) 


Combining  Eq.(43),  Eq.(45)  and  Eq.(48)  and  further  Eq.(43),  Eq.(44)  and  Eq.(48)  we 
obtain: 


axb3  -a3bx  =  0  , 


(49) 


axb2  +  a2bx  =  0  . 


(50) 
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Note  that  relation  Eq.(47)  can  be  obtained  by  combining  Eq.(44),  Eq.(45)  and  Eq.(48)  as 

well. 

For  further  practical  usefulness  the  relations  Eq.(47),  Eq.(49),  Eq.(50)  can  be  rewritten  in 
the  form: 

axb2  -  a2bx  =  0  , 

( axb2  +fl2Z)1)sin2(a-(p)+(fl2fl3  +  Z>263)cos2(a-cp)  =  0,  (51) 

(axb2  +  a2bx  )cos  2(a  -  (p )-(a2a3  +  b2b3  )sin  2(a  -  (p)  =  0  . 

Finally,  substitution  Eq.(46)  in  Eq.(51)  gives  the  relations  for  determining  the  values  of 
anisotropy  parameters  for  medium  having  the  orthogonal  eigenpolarizations: 

(1  -  P)2  R cos2  y  -  (1  -  P2 )(1  -  R2 )  sin ^ cos ^ sin  2{a  -  y)  -  (1  +  P)1  R sin2  §  =  0 ,  (52) 

(1  -  P)[(l  +  R 2 )  cos  2(a  -  y  -  </>)  -  (1  -  R 2 )  cos  2{a  -  y)\  =  0 ,  (53) 

2(l  +  P2  V?  sin  —  cos  —  + 
v  ’  2  2 

+  ^(l-/>2)sin2  y[(l-7?2)sin2(a-y)  +  (l  +  i?2)sin2(«-y-^»)]-  .  (54) 

-  (1  -  P2) cos 2  y  [(1  -  R 2 )  sin  2 (a  -  y)  -  (1  +  R 2 )  sin  2 (a  -  y  -  <p)\  =  0 


From  Eq.(52)  it  is  follows=> 


A  =  -2  arctg 


1  P- 1 


2  R  P  + 1 


(1  - R2 ) sin 2 (a  -/)- -^/sin2  2(a-y)(l-7?2)2  +47?" 


(p  =  a  -  y  +  —  arccos|  cos2(o:-y) 


R2  -1 
R2  +1 


n 

~2 


(55) 


(56) 


Eq.(52)  and  Eq.(54)  can  be  represented  as 

-  2PR  +  (l  +  P1  )r  cos  A  +  ^  (l  -  P2  )(l  -  R2  )sin  A  sin  2{a  -  y)  =  0 , 

(l  +  P2  )r  sin  A  +  -^  (l  -  P 2  )[(l  +  R2  )sin  2  (a  -  y  -  (p)-{[-  R2  )sin  2(a  -  y)  cos  a]  =  0 . 
Making  the  following  designations 

A  =  ^(l-P2\\-R2)sm2(cc-y);  B  =  (\  +  P2)r ; 

C  =  -2  PR;  D  =  -^(\-P2\\  +  R2)  sin  2  {a-y-  cp), 

then  Eqs.(57)  take  the  forms 


(57) 


(58) 
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^4  sin  A  +  i?  cos  A  +  C  =  0, 
^4  cos  A  -  5sin  A  +  D  =  0. 


(59) 


For  Eq.(59)  to  be  satisfied  the  following  equality  is  required 

A2  +B2  =C2  +D2 . 

Taking  to  account  Eqs.(59)  and  (60)  we  have: 

(C  —  B)tan2  —  +  2  A  tan  —  +  ( C  +  B )  =  0 , 

(D  -  A) tan2  y  -  2.5  tan  y  +  (Z)  +  A)  =  0. 
Solution  of  the  system  Eqs.(61)  is 

A  -zf±V^42  -C2  +52  -A±4d2  -52  +  52 

tg—  = - = - 

2 1,2  C-B  C-B 


A  B  +  ^B^CC^A2  BC^B2  ^B^+C2 

tg—  = - = - 

2 1,2  D- A  D- A 

From  Eq.(60)  it  follows: 

A+D _B+C 
C-5  _  Z)-zf 


-A  +  D 
C-B 


l  C-B 
B  +  C 
D-A 
B-C  ' 

. D-A 


A  +  D  ^  B-C 
C-5  *  D^A' 

Evidently,  we  need  to  choose  in  Eqs.(62)  and  (63)  the  same  solutions,  then 

A  A  +  D  A  -A  +  4a~2  -C2  +52 

tg  — = - =  tg  — = - . 

2  B-C  2  C-B 


(60) 


(61) 


(62) 


(63) 


(64) 


(65) 


As  it  can  be  seen  Eq.(65)  is  equal  to  Eq.(55).  Thus,  we  have  2  equations  and  6  unknown 
parameters  to  be  determined.  Therefore,  4  out  of  6  parameters  we  should  fix  and,  then,  other  2 
can  be  determined  from  Eq.(52)  and  Eq.(53).  It  turns  convenient  to  determine  the  parameters  A 
and  (p  (see  Eqs.  (55)  and  (56)),  while  other  4  parameters  are  fixed  in  a  range  of  definition. 

For  depolarizing  Mueller  matrix  interpretation  the  method  of  coherency  matrix  introduced 
by  Cloude  [35,36]  has  widely  been  used.  According  to  this  method  an  arbitrary  Mueller  matrix 
M  including  depolarization  and  errors  effect  are  represented  as  a  sum  of  four  nondepolarizing 
(deterministic)  Muller  matrices  M^4 : 


k=  1 


(66) 


In  Eq.(66)  ///f  -  are  eigenvalues  of  the  Cloude ’s  coherency  matrix  playing  the  role  of 
weighting  multipliers.  So,  in  compliance  with  Eq.(66)  the  anisotropic  properties  of  object  are 
presented  by  simultaneous  parallel  independent  effects  of  four  deterministic  anisotropic  parts. 

If  in  Eq.(66)  only  one  eigenvalue  /4  is  non-zero,  the  corresponding  Mueller  matrix 

16 
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represents  nondepolarizing  (homogeneous  anisotropic)  medium.  If  all  four  eigenvalues  of 
coherence  matrix  are  non- zero  the  dominant  type  of  deterministic  effect  of  object  on  input 
polarization  is  described  by  Mueller  matrix  associated  with  the  maximum  eigenvalue  in  sum 
Eq.(66).  The  case  of  negative  eigenvalues  indicates  that  the  Muller  matrix  describes  a  physically 
impossible  transformation.  This  may  occur  as  a  result  of  measurement  errors.  If  eigenvalues  of 
coherency  matrix  are  comparable,  then  the  studied  medium  depolarizes  input  radiation 
completely.  In  this  case  dominant  type  of  deterministic  behaviour  of  studied  medium  can  not  be 
determined. 

Our  previous  studies  have  shown  that  direct  using  of  the  generalized  equivalence  theorem 
for  the  analysis  of  “raw”  experimental  Mueller  matrices  can  lead  to  significant  errors  in  the 
results  analysis.  This  means  that  before  interpreting  initial  experimental  Mueller  matrix  should 
be  fdtered  out  to  exclude  negative  eigenvalues  of  coherency  matrix. 

The  elements  of  Cloude’s  coherency  and  Mueller  matrix  are  related  as  follows  [35,36]: 


c 


c 


c 

c 


11 


12 


13 


14 


‘■'31 


'32 


'33 


'34 


4 

4 

4 

J_ 

4 

J_ 

4 

4 

J_ 

4 

J_ 

4 


(mn+m22+m33+m44),  c21=- 
(ml2  +  m21  -  im34  +  im4i ) ,  c22  = 

(ml3+m3l+im24-im42),  c23  = 

(mX4  - im23  +  im32  +m4l),  c24  = 

(ml3  +  m3l-  im24  +  im42 ) ,  c4l  = 
( -im14  +  m23  +  m32  +  im4X  ) ,  c42 

(mn-m22+m33-m44),  c43=- 

(imn  - imlx  +  m34  +m43),  c44  = 


-(ml2+m2l+im34-im43), 

=  ^(mn+m22-m33-m44), 
--^(imH+m23+m32-im4l), 

=  ^ -[-iml3  +  im3X  +  m24  +  m42 ) , 

--^(ml4+im23-im32+m4l), 

=  ^ ( imx 3  - imiX  +  m24  +  m42 ) , 
^(-/m12+/m21+m34  +  m43), 
^(mn-m22-m33+m44). 


(67) 


Coherency  matrix  C  has  always  four  real  eigenvalues,  since  it  is  Hermitian,  i.e.,  ct  .  =  c*  j , 

where  *  denotes  Hermitian  conjugation.  Four  eigenvectors  'P1-4  of  the  matrix  C  are  directly 
related  to  the  Jones  matrices  T1-4  as  follow: 


f(k)  \T/(^)  ,  \T/(^)  Ak)  \T/ (k)  _ 

^11  T1  "r  T  2  5  L12  T3  T4  ’ 

f(k)  \T/ (k)  ,  j\jf(k)  Jk)  \T/ (k)  _  \T/ (k)  . 

L2l  T  3  '  L  T  4  5  L22  T  1  T  2  ’ 


k  =  1-4; 


(68) 


where  k  -  number  of  eigenvectors  of  coherency  matrix  C . 

The  initial  depolarizing  Mueller  matrix  is  represented  by  the  sum  of  four  deterministic 
Muller  matrices  M'j4  corresponding  to  the  Jones  matrices  Eq.(68). 

Eigenvalues  of  the  coherency  matrix  C  are  used  to  characterise  the  depolarization 
properties  of  the  studied  object  by  single-value  depolarization  metric  called  entropy  H. 
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where  N  =  1,4  is  chosen  such  that  the  condition  0  <  H  <  1  is  satisfied.  In  the  case  when  all 
eigenvalues  of  the  coherency  matrix  are  positive  one  should  set  N  =  4 .  In  other  cases,  N  is  the 
number  of  positive  eigenvalues  of  the  coherency  matrix  C. 

Entropy  H  defined  in  Eq.(69)  characterizes  the  degree  of  anisotropic  disorder  yielding 
depolarization  of  incident  radiation.  If  H  =  0 ,  i.e.,  only  one  eigenvalue  of  coherence  matrix  is 
nonzero,  the  studied  medium  described  by  corresponding  deterministic  Mueller  matrix  is 
homogeneous  and  doesn't  depolarize  input  radiation.  If  H  =  1 ,  the  studied  medium  is  ideal 
depolarizer.  Intermediate  values  of  H  correspond  to  anisotropic  depolarizing  medium  that  in 
addition  to  the  depolarization,  change  the  polarization  of  input  radiation.  As  it  was  noted  above, 
the  author  [35,36]  suggested  that  a  medium  leading  to  depolarization  can  be  characterized  by  the 
some  dominant  effective  anisotropy  in  case  of  H  <  0.5 .  If  in  Eq.(66)  all  four  components  have 
comparable  weighting  multipliers,  all  of  them  contribute  equally  to  the  change  of  polarization  of 
radiation. 

2.4  Development  of  optimized  polarimeter 

In  given  section  we  describe  an  assembled  polarimeter  that  is  optimized  for  accurate 
measurement  of  all  elements  of  Mueller  matrix  of  media.  As  we  noted  early  it  would  be 
extremely  desirable  that  this  polarimeter  could  dynamically  adjust  to  exact  experimental 
conditions  and  have  no  moving  units.  This  will  provide  a  minimal  measurement  time  with 
maximum  accuracy  and,  hence,  be  a  promising  basis  for  new  type  of  imaging  Mueller 
polarimetry.  One  of  possible  ways  to  realize  dynamically  adjusted  polarimeter  which  we 
implemented  here  is  using 

LC  retarders  as  polarization  transducers,  which  retardation  is  electrically  variable  -  Liquid 
Crystal  Variable  Retarders  (LCVR).  As  it  was  shown  in  [16]  a  minimal  number  of  LC 
transducers  to  measure  complete  Mueller  (all  of  sixteen  matrix  elements)  matrix  is  four. 

Mentioned  LC  transducers  (Fig.l.)  can  replace  an  entire  series  of  polymer  and  ordinary 
crystalline  retarders.  With  new  Swift  LC  Meadowlark  technology  (www.meadowlark.com),  the 
switching  speed  is  symmetric  and  approximately  150  microseconds. 


Fig.l.  Meadowlark’s  Liquid  Crystal  Variable  Retarder  (LCVR). 
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The  Four  Channel  Digital  Interface  D3050  (Fig.2)  designed  for  computer  control  of  up  to 
four  Meadowlark  Optics  nematic  liquid  crystal  devices  at  the  same  time.  Package  allows  the 
amplitude  of  the  2  kHz  square  wave  output  to  be  driven  either  by  an  external  DC  analog  signal 
supplied  to  a  front  panel  connector  or  specific  CellDRIVE  (Fig.3)  generated  waveforms 
including  sinusoidal,  square,  triangle,  sawtooth  and  transient  nematic  effect  waveforms. 
Additional  functions  include  the  capability  to  output  a  sync  pulse  on  a  front  panel  connector  at 
desired  points  in  the  CellDRIVE  generated  waveforms  and  the  ability  to  save/restore  all 
CellDRIVE  settings  to/ffom  a  file. 


Fig.2.  Meadowlark’s  Four  Channel  Digital  Interface. 


- CELLDRIVE  3100  ,+Wr 
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1  ituml  Eaii™ 

EEE  1329  im 


pn1.'  =guinFiiiii=gj^^M  i. 


Sync  Gul  Solaria 

Gui  Fhasr  (draj  r-oui 


uiiiniirir 


TimpUpto 


Fig.3.  CellDRIVE  3100  advanced  software 


The  tests  of  purchased  LC  transducers  (for  the  results  see  Fig.4)  was  realized  on 
automatized  Mueller-polarimeter  which  operates  with  rotatable  crystal  plates  at  input  and  output 
channels  (Fig.  5). 

Fig.4  shows  nonzero  anisotropies  of  LC  transducers  during  applied  voltage  changing.  It 
can  be  seen  that  besides  main  type  of  anisotropy  (linear  birefringece  -  A)  the  retarder’s  cell  have 
also  dependence  of  difference  in  transmittance  for  eigenpolarizations  (linear  dichroism)  versus 
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applied  voltage.  This  is  very  important  fact  which  must  to  be  taken  in  to  account  during 
polarimeters  designing. 
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Fig.4.  Dependences  of  LC  transducers  parameters  on  value  of  applied  voltage.  Analysis  is 
performed  in  conformity  with  theorem  Eq.(10). 


Fig.5.  Scheme  of  Mueller  matrix  polarimeter  used  for  LC  transducers  tests. 

1.  He-Ne  laser; 

2.  polarizer; 

3.  rotatable  crystal  phase  plate  which  orientation  during  measurement  periodically 
consequently  takes  the  discrete  values  (a/); 

4.  LC  transducers  tested; 
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5.  continuous  rotatable  crystal  phase  plate; 

6.  analyzer; 

7.  photodetector; 

8.  PC; 

9.  digital  interface  D3050. 

Figure  6  shows  the  scheme  of  Mueller-polarimeter.  Receiving  channel  (PSA)  is  a  complete 
Stokes-polarimeter  using  two  LC  transducers  (LCVR3  and  LCVR4),  described  above.  In  the 
probing  channel  (PSG)  consists  of  two  LC  transducers  (LCVR1  and  LCVR2)  as  well. 
Orientations  of  all  units  in  polarimeter  scheme  are  indicated  in  Fig.6.  Mathematics  of  polarimeter 
operation  is  presented  below. 
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where  M  is  a  Mueller  matrix  of  object,  G  -  characteristic  matrix  of  Mueller-polarimeter,  r‘  s', 
are  /-th  Stokes  parameters  for  /-th  polarization  of  input  and  transmitted  light  respectively. 


Stokes-polarimeter 
(Polarization  state  analyzer  (PSA)) 


Fig.6.  Mueller-polarimeter  based  on  four  LC  transducers. 

I .  He-Ne  laser; 

2,8.  polarizer  (analyzer)  with  orientation  6P  ( 0A ); 

3, 4, 6, 7.  LC  transducers  (LCVR1-4)  which  orientations  of  fast  axis’s  are  aw 
and  introduced  additional  phase  shifts  are  A1.4  (indexes  i,j  denotes 
changed  phase  shifts); 

5.  test  sample; 

9.  photodetector; 

10.  PC; 

II.  digital  interface  D3050. 

Ii  -  intensity  registered  by  detector; 


Distribution  A:  Approved  for  public  release;  distribution  is  unlimited. 


21 


R,  S  -  Stokes  vectors  of  input  and  transmitted  light  correspondingly; 

M  -  Mueller  matrix  of  measured  samples. 

According  to  Eq.(70),  measurement  of  the  Mueller  matrix  by  polarimeter  Fig.6  is  occurred 
by  the  means  of  sequential  irradiation  of  sample  by  electromagnetic  radiation  with  four  different 
polarizations. 

Through  the  fact  that  in  Eq.(70)  the  Stokes  parameters  of  output  radiation  are  the  measured 
parameters  it  needs  to  consider  a  PSA  of  polarimeter,  Fig.6,  in  more  details.  PSA  is  a  Stokes- 
polarimeter  based  on  two  LC  transducers  (see  scheme  Fig.7). 


Stokes-polarimeter 
(Polarization  state  analyzer  (PSA)) 


LCVR-4  lcvr.3 


Fig.7.  Stokes-polarimeter  based  on  two  LC  transducers. 
LCVR-3,4  -  LC  transducers; 

A-  analyzer; 

D  -  detector; 

S  -  Stokes  vector  of  radiation  to  be  analyzed. 


Expression  for  intensity  (i.e.,  first  parameter  of  the  Stokes  vector)  after  analyzer  A,  Fig.7, 
can  be  written  as 

I  =  (M,  (0°)  •  Mlcvr_4  (A4  ,«/)•  Mlcvr_3  (A3 ,  a3 )  •  s)( .  (71) 


where  M  A{0) ,  MLCVR_3  4  (A,  a)  are  the  Mueller  matrices  of  analyzer  with  orientation  6  and  LC 

transducers  with  phase  shift  8  and  orientation  a  . 

Due  to  the  orientation  a  of  LC  transducer  can  not  be  changed  electronically  for  measuring 
of  all  four  Stokes  parameters  it  is  necessary  to  change  the  values  of  birefringence  of  LC 
transducers  four  times  A{'“4)  at  least.  Then,  we  can  set  a  system  of  equations  relative  to 
measured  Stokes  parameters: 


A 


i  = 


A 

A 

A 


(M4(0o)-M  lcvr_4  (A« ,  «4 )  •  M  lcvr_3  (A«  ,a3)-s\ 
(M  4  (0°)  •  M  LCVR_4  (Af ,  a4 )  •  M  LCVR_3  (A<2)  ,a3)-S\ 
(M4  (0°)  •  M  LCVR_4  (A<3) ,  a4 )  •  M  LCVR_3  (A<3) ,  a3 )  •  S  )( 
(M  4  (0°)  •  M  LCVR_4  (A<4) ,  a4 )  •  M  LCVR_3  (A*4)  ,a3)-S\ 


I  =  B(A3\a3,A4\a4)  ■  S  .  (72) 


It  is  known  [31]  that  the  system  Eq.(72)  can  be  solved  with  minimal  error  when  its 
condition  number  VB  is  minimal: 

VB  =||  B(A«,a3,A«,a4)  ||  ||  B( A« , cr3 , A? , cr4 )"‘  ||^  min  =>  ->  ^Smin.  (73) 

where  ||  ||  denotes  Euclidian  (metric)  norm  and 
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5  S  = 


(74) 


AS 


S  =  S0  +  AS . 


In  Eq.(74)  S0  is  exact  (table)  Stokes  vector;  AS  denotes  vector  of  errors  of  corresponding 
Stokes  parameters. 

Eq.(73)  results  from  the  next  relation  [31] 


SS< 


VB  -(<SB  +  <3I) 

1  -  V^SB  ’ 


(75) 


where  SB,  SI  -  are  errors  of  determination  of  matrix  B  (called  characteristic  matrix  of  the  system 
of  linear  equations  or  instrumental  matrix  of  Stokes-polarimeter)  and  of  vector  I  of  intensities 
measured  by  detector  D,  see  Fig.  7. 

We  have  found  that  the  minimum  value  of  condition  number  is  Vmia  =  4.5  and  it  is  reached 
when  angle  between  orientations  of  LC  transducers  is 

|a4  -a3|  =  22.5°  (76) 


and  its  phase  shifts  are 

A<'-4)  =  (3.2°,  162.5°,  32.6°,  167.0°); 
A<1_4)  =  (131.1°,  53.4°,  1.3°,  159.7°  ); ' 


(76') 


Fig.  8.  Dependence  of  condition  number  on  orientations  of  LC  transducers  a3 ,  a4. 


After  determination  of  actual  range  of  phase  shift  changing  for  LC  transducers 
(A3,4e[30o,180°]),  we  corrected  the  conditions  for  optimal  determination  of  Stokes  parameters  by 
polarimeter  Fig.7.  As  a  result  the  optimal  set  of  values  of  phase  shifts  of  LC  transducers  in 
comparison  with  previous  estimations  has  changed.  New  (corrected)  set  is 
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(77) 


A<'-4)=(40,  40,  180,  180°); 
A^-4)=(40,  138.52,  173.26,  59.90  ) 


At  the  same  time,  the  optimum  orientation  of  LC  transducers  remains  the  same  as  previous 
one  Eq.(76) 

When  the  optimum  set  of  parameters  of  transducers  was  established,  we  analyze  how  the 
imperfections  of  these  transducers  and  other  elements  of  polarimeter  affect  an  accuracy  of  each 
of  the  four  Stokes  parameters.  To  estimate  errors  for  each  Stokes  parameters  we  used  the  known 
relation  for  the  errors  of  indirect  measurements  [31]: 


D-P(ki,y2 


yN)  = 


dF(y°K) 

dyn 


y 

) 


(78) 


where  F(yl,y2,...,yN)  designates  some  function  of  measured  variables  yn;  y°n  -  mathematical 
expectation  (averaged  value). 

According  to  Eqs.(72)  and  (78)  the  expression  for  the  estimations  of  Stokes  parameter’s 
errors  take  the  form: 
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(79) 


To  obtain  quantitative  estimations  for  expression  (79)  we  set  the  next  values  of  errors: 
a«=0.2°,  aA=0.5°  and  a/=0.1%  (taking  the  intensity  as  /=  1).  As  a  result,  basing  on  Eq.(79),  we 
obtain  the  following  dependences  of  errors  for  each  of  Stokes  parameters  aS,.  on  azimuth  [3  and 

ellipticity  angle  s  of  analysed  polarization  ellipse  (Fig.9). 
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Fig.9.  Dependences  of  Stokes  parameter’s  errors  on  polarization  of  analysed  radiation  (/?  - 
azimuth  of  polarization  ellipse,  s  -  ellipticity  angle). 

Mean  AS  and  peak  z\S/<;,  value  of  error  for  each  Stokes  parameter  is: 

AS  =  [0.005  0.007  0.010  0.008]  ^ 

AS^  =  [0.001  0.001  0.006  0.003] 

As  it  is  for  Stokes-polarimeter  the  optimal  polarizations  R<1_4)  of  input  radiation  are 
generated  for  the  next  pairs  of  values  of  phase  shifts  of  LC  transducers  1  and  2: 


A[’-4)  =(40,40,180,180°); 

A<'-4)  =  (40,138.52,173.26,59.90  ); 


This  result  was  obtained  by  minimisation  of  condition  number  of  the  characteristic  matrix 
G  Eq.(70).  The  expressions  for  the  Stokes  parameters  of  input  radiation  are: 

rn 

R“’  =  ( Af, a, ).MIOTM(Af >,).  ’  .  (82) 

v°. 

In  accordance  with  scheme  Fig.6  we  have  assembled  the  dummy  of  Mueller-polarimeter 
and  have  developed  software  modules  to  control  of  its  operating.  Improved  software  interface  is 
shown  below  in  Fig.  10. 
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Panels  for  representation  of  the  results  of  Stokes  Panel  for  settings  the  parameters  of  rotating 

parameters  and  Mueller  matrices  measurements  phase  plate  and  ADC  parameters 


PI 


P2 


-I6M  x| 


Panel  with 
results  of 
S  application  of 
Cloude 

decomposition 
and  Generalized 
equivalence 
theorem 
P6 


Panel  for  driving  of  mueller-polarimeter  with  4  LC  transducers 

P4 


Panel  for  settings  the  LC  transducers  parameters 
P5 


Fig.  10.  Developed  software  for  Mueller  matrix  polarimeter  operation  and  control  which  realizes 
the  time-sequential  strategy  with  using  four  LC  transducers.  Modules  with  Cloude’s 
decomposition  and  generalized  equivalence  theorem  are  implemented  as  well  (see  panel  P6). 


Measurement  procedure  with  developed  polarimeter  consists  of  the  next  sequence  of 
operations: 

1.  First,  we  need  to  define  polarizations  of  probing  radiation.  To  do  this  we  need  to  remove 
all  objects  from  probing  channel  of  polarimeter  (measurement  mode  “without  object”)  and  press 
the  button  “Calibrating  vs  LC”  of  developed  software  which  placed  on  driving  panel  P4  (see 
Fig.  10).  Then,  one  has  to  follow  to  options  of  corresponding  dialog  menu.  As  a  result  input 
Stokes  vectors  are  formed  and  represented  on  the  right  of  panel  PI. 

2.  Second,  to  measure  the  Mueller  matrix  we  insert  an  investigated  object  on  the  way  of 
laser  beam  after  PSG  (see  Fig.6)  and  press  the  button  “Matrix  vs  LC”  and  follow  to  the  options 
of  corresponding  dialog  menu  again.  Measured  matrix  will  be  figured  on  panel  PI  as  M: 
(nonnormalized  one  is  on  the  top  and  normalized  by  mu  on  the  bottom). 

3.  If  we  want  to  determine  the  parameters  of  polarization  (<9  -  azimuth,  s  -  ellipticity,  I  - 
intensity,  P  -  polarization  degree,  S  -  Stokes  parameters)  which  measured  by  PSA  (Fig.6)  we 
have  just  press  the  button  “Stokes  vs  LC”  and  wait  until  required  information  will  be  figured  on 
the  corresponding  panel  PI  (see  Fig.  10). 

3.  Results  and  Discussion 

Let  verify  the  derived  conditions  of  eigenpolarizations  orthogonality  for  some  simple 
cases.  We  will  analyze  the  orthogonality  conditions  in  terms  of  the  Mueller  matrix  elements 
Eq.(31),  Eq.(42)  (condition  Eq.(28)  is  omitted  because  it  is  a  particular  case  of  Eq.(42)). 

As  it  was  pointed  out  previously,  the  media  characterized  by  only  one  of  the  basic  types  of 
anisotropy  possess  always  orthogonal  eigenpolarizations.  Mueller  matrices  of  these  types  of 
anisotropy  are  given  in  Eq.(12).  As  can  be  directly  seen,  for  all  of  them  the  condition  Eq.(41)  is 
always  satisfied.  Conditions  Eq.(31)  in  these  cases  of  anisotropy  become  an  identity  or  have  the 
uncertainty  of  the  type  0/0. 
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Multiple  numerical  verifications  show  that  if  the  model  Eq.(10)  consists  all  anisotropy 
parameters,  conditions  Eq.(31)  and  Eq.(42)  are  satisfied  if  values  of  circular  and  linear 
anisotropy  satisfy  Eq.(55),  Eq.(56).  For  example: 


P  =  0.2;  y  =  25°;  R  =  0.7;  a  =  84°;  =>  <p  =  9.38°;  A  =  -51.83°.  (83) 


T  = 


rp  CP  rp  IP  rp  CA  rp  LA 


0.803  +  0.61/  0.457  +  0.098/ 

0.188  +  0.428/  0.261  +  0.168/ 


(84) 


It  is  easy  to  see  that  eigenpolarizations  X\,i  °f  the  Jones  matrix  Eq.(84)  are  orthogonal. 
Indeed: 


Xx  =  0.446  +  0.228/;  =  -1 .777  -  0.908/; 

H5  't 

Z1Z2  =-!• 


Taking  to  account  Eq.(83),  in  terms  of  Mueller  matrix  we  obtain: 


M  =  MCPM/PMC4MZ“4 


1.000 

0.593 

0.707 

0.361 


0.593 

0.436 

0.356 

0.229 


0.707 

0.394 

0.568 

0.197 


0.361 

0.155 

0.260 

0.238 


Eigenpolarizations  of  the  Mueller  matrix  M  Eq.(86) 


"1.000" 

"  1.000  " 

T 

0.599 

-0.599 

0 

s,  = 

;  S2  = 

;  S,  +S2  = 

I 

0.713 

-0.713 

0 

0.365 

-0.365 

0 

IlpH  3tom  Eq.(31)  =>  0.034=0.034=0.034 


(85) 


(86) 


(87) 


In  order  to  check  the  orthogonality  conditions  in  terms  of  anisotropy  parameters  Eq.(52)- 
Eq.(54)  let  us  consider  some  special  cases. 

(i)  In  case  of  absence  of  amplitude  anisotropy  ( P  =  1  and  R  =  0 )  we  have: 


0  =  0 

0  =  0.  (88) 
0  =  0 


The  Jones  matrix  for  the  case  of  P  =  1  and  R  =  0  is  a  unitary  [2]  and,  indeed,  characterized 
by  orthogonal  eigenpolarizations. 

(ii)  In  the  case  of  absence  of  phase  anisotropy  ( A  =  0  and  <p  =  0 ),  Eq.(28)  takes  the  form: 

(l-P)2P  =  0 

2(1  -  P)i?2  cos2y  =  0  (89) 

(1-P2)P2  sin 2y  =  0 
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This  proves  again  that  simultaneous  presence  of  both  type  of  amplitude  anisotropy  results 
in  non-orthogonality  of  eigenpolarizations  [1], 

It  is  important  to  note  that  relations  Eq.(52)-Eq.(56)  can  be  considered  as  a  set  of  rules  for 
synthesis  of  the  polarization  elements  with  orthogonal  eigenpolarizations  as  well. 

Results  obtained  in  section  2.1  show  that  Mueller  matrix  of  the  medium  with  orthogonal 
eigenpolarizations  has  characteristic  symmetry,  at  least  such  as  in  Eq.(42).  At  the  same  time  we 
show  that  such  class  of  media  can  generally  contain  all  types  of  anisotropy. 

Next  important  result  is  that  the  medium  with  orthogonal  eigenpolarizations  have  no  more 
than  4  degrees  of  freedom.  Thus,  we  can  claim  that  the  knowledge  of  incomplete  Mueller 
matrices  Eq.(90)  are  sufficient  for  complete  description  of  anisotropy  of  the  medium  with 
orthogonal  eigenpolarizations: 


mn 

* 

* 

* 

mn 

mn 

ml3 

1 

3" 

m2l 

m22 

m23 

m24 

* 
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m33 

m34 
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3 

4^ 

m42 

m43 

3 

4^ 

4^ 

1 _ 

* 

m42 

m43 

3 

4^ 

4^ 

1 _ 

(90) 


Furthermore,  taking  into  account  the  relations  [34]  the  number  of  matrix  elements  that 
need  to  be  measured  may  be  further  decreased. 

Fig.4  shows  that  LC  transducers  in  addition  to  linear  birefringence  has  an  evident  linear 
amplitude  anisotropy  depending  on  applied  voltage  as  well.  This  needs  to  be  taken  into  account 
during  polarimeter  adjustments.  In  particular  it  is  reasonable  to  operate  with  applied  voltage 
more  then  3V. 

It  was  shown  in  [14]  that  minimum  value  for  condition  number  of  instrumental  matrix  with 
dimension  4x4  for  time-sequential  Stokes-polarimeter  with  rotating  phase  plate  which  can 
principally  be  obtained  is  Vmin  =  4.47 .  For  the  case  of  Stokes-polarimeter  with  two  fixed  LC 
transducers  we  obtain  value  Vmin  =  4.50  that  is  not  far  from  principally  obtainable  minimum. 
This  allows  providing  minimal  errors  in  measurement  of  Stokes  parameters.  Note,  these  results 
concerns  so-called  integral  errors  of  Stokes  vector  measurement  Eqs.(74)  and  Eq.(75). 

Unfortunately,  in  practice  an  actual  range  of  change  of  phase  shift  for  LC  transducers  lies 
in  region  A  e  [30°,  180°]  (smooth  regions  on  dependencies  presented  in  Fig.4).  So,  we  have  found 
new  set  of  A  1,2,3, 4  Eq.(77).  This  results  in  some  increasing  of  condition  number  Emm  (from  4.5  to 
4.6).  However,  under  equal  other  conditions  this  increases  the  limits  of  relative  error  dS  just  up 
to  1 .02  times. 

As  it  can  be  seen  from  Fig.9,  the  error  of  Stokes  parameters  determination  by  polarimeter 
Fig.7  depends  on  the  polarization  of  input  radiation.  This  dependence  is  symmetrical  relative  to 
sign  of  ellipticity  of  polarization  ellipse  for  parameters  y  and  s3,  and  is  asymmetric  for 

parameters  s2  and  s4.  All  dependences  in  Fig.9,  except  for  parameter  s3,  are  asymmetric 

relative  to  zero  value  of  azimuth  J3.  From  Eq.(80)  it  can  also  be  concluded  that  in  the  scheme  of 
polarimeter,  which  is  considered,  the  Stokes  parameters  are  generally  measured  with  different 
accuracy.  Meanwhile  parameter  y  is  most  precisely  defined,  s2  a  little  worse,  than  S4and 

parameter  s3  defined  with  lowest  precision.  For  the  parameter  s3  the  magnitude  of  error  is  also 

highest.  The  positions  of  extremum  of  dependences  Fig.9  are  determined  by  LC  transducers'  axes 
orientation  (  ax  =  22.5° ,  a2  =  45°  ). 

Thus,  in  considered  polarimeter  scheme  for  arbitrary  polarization  of  radiation,  intensity 
will  be  determined  with  maximum  accuracy.  Polarization  parameter  -  large  axis  of  the 
polarization  ellipse  that  oriented  at  0°  and  90°,  will  also  be  determined  with  the  greatest  accuracy 
for  arbitrary  ellipticity.  Accuracy  of  analysis  of  the  polarization  ellipse  orientation  close  to  45°  is 
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the  worst,  and  strongly  depends  on  ellipticity  of  ellipse.  Ellipticity  angle  of  polarization  ellipse  in 
described  scheme  of  polarimeter,  taking  into  account  comparability  of  values  Us4  and  U s  4  is 

determined  with  greater  error  relative  to  azimuth  determination  and  depends  on  the  orientation  of 
the  polarization  ellipse. 

The  algorithms  of  Cloude  decomposition  and  decomposition  in  accordance  with 
generalised  equivalence  theorem  described  in  section  2.2  are  implemented  in  appropriate 
software  module  Fig.  10.  It  provides  additional  information  about  studied  objects.  In  particular,  in 
the  case  shown  in  Fig.  10,  it  can  be  seen  that  measurement  of  Mueller  matrix  was  carried  out  in 
the  mode  “without  object.”  In  this  case  the  measured  Mueller  matrix  is  unitary  diagonal  within 
error  range  8M  «1.9%.  This  do  result  in  relatively  large  entropy  H=  0.11.  Values  of  anisotropy 
parameters  indicate  that  measured  object  (empty  space)  does  not  characterized  by  any  type  of 
anisotropy  just  as  it  had  expected.  Next  figure  demonstrates  the  part  of  interface  Fig.  10  when 
prism  polarizer  is  measured: 


Results  of  analisys= 
Experimental  matrix 


0,274502 
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Fig.  1 1 .  Part  of  interface  Fig.  10  demonstrating  the  case  when  measured  object  is  a  prism 

polarizer. 

Fig.l  1  represents  the  Mueller  matrix  with  structure  typical  for  linear  amplitude  anisotropy: 
entropy  H  is  about  0.08;  main  anisotropy  type  (linear  amplitude  anisotropy  with  value  “.P  -0, 
and  azimuth  “Teta”  =  -37°. 

Analysis  of  the  Muller  matrix  of  structure  Eq.(90)  has  shown  that  its  elements  can  not  be 
determined  independently  by  means  of  time-sequential  approach  described  above,  see  Eq.(70). 
Since  the  first  row  of  the  matrix  is  responsible  for  converting  the  intensity  of  the  radiation  by 
medium,  they  will  be  involved  in  the  construction  of  relations  Eq.(70)  [37,38], 

From  the  other  hand,  for  determining  of  elements  of  the  Mueller  matrix  there  is  no  need  to 
use  the  Stokes  polarimeter  in  receiving  channel  of  polarimeter.  System  of  simultaneous 
equations  relative  to  13  matrix  elements  of  the  form  Eq.(90)  can  be  composed  basing  on 
relations  Eq.(72)  directly: 


>  A(t) 

?  ^4  ? 

a3)  M13 

•S(A<4),A  [k),a 

II 

<N 

+  r2Sl 

)mn 

+  0^3  + 

r3sx)mn  +  (rrs4 

+  r4sx)ml4  + 

+r2s2 

m22 

+  r2s3m23 

+  r2s4m24  + 

+r3s2 

m  3  2 

+  r3s3m33 

+  r3s4m34  + 

+r4s. 

2^42 

+  r4s3m43 

+  r4s4m44 , 

(91) 


where  k=l-\3,  It  -  measured  intensity  in  the  case  of  k- th  set  of  phase  shifts  of  LC  transducers 
and 
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R(^,Ai*>,A<*»,a4,a3)  =  [M,(^)-MIOT_4(A«>,a))-MIOT_3(Af,aJ)]  , 

1 


S(^ ,  A<‘> ,  A<‘> ,a2,ax)  =  MLCVR_2  (A^ ,  a2 )  •  MLCVR_X  (A}*' ,  a, ) 


(k) 


iW 


cos(2#p) 

sin(2#p) 

0 


(92) 


<1>  -  denotes  first  row  of  resulting  matrix  and  MJ3  denotes  the  following  Mueller  matrix: 
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From  Eq.(91)-  Eq.(93)  we  can  form  the  system  of  equations: 


(93) 


Mi3;  it  =  1-13.  (94) 


B13x13 

where  I  -  vector  of  intensities,  Bi3xi3  -  characteristic  matrix  of  dimension  13x13;  M13  -  the 
Mueller  vector  of  the  form: 
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From  Eq.(94)  we  have  for  mij  : 


M.3=(B13i1j)-'i.  (96) 

As  it  was  in  the  case  of  Stokes-polarimeter,  see  Eqs.(72)-(75),  to  minimize  an  error  of 
Mueller  matrix  elements  determination  in  Eq.(94)  we  have  minimize  the  condition  number 
VR  of  characteristic  matrix  Bi3xi3.  From  Eq.(94)  it  can  be  seen  that  V„  is  a  function  of  10 

parameters  (0A,0p,al_4,Al_4).  It  is  obvious  that  we  can  set  0P  =  0  and  other  azimuths  are 
determined  relative  to  0P .  Then  ^  was  numerically  minimized  taking  into  account  that  actual 

range  of  phase  shift  changing  for  LC  transducers  is  Ai_4e[30°,180°].  As  a  result  the  minimum 
value  V„  .  w  20.23  under  next  conditions: 

^13*13  >mm 


-ao-13)- 

'42.4°  91.5°  122.8°  40.0°  135.4°  180°  120.4°  40.9°  180°  177.8°  73.1°  40.5°  84.5°  ~ 

a(2'-13) 

54.7°  179.8°  174.3°  177.4°  179.9°  54.8°  145.0°  178.5°  40.1°  41.0°  44.6°  173.8°  65.6° 

ao-13) 

95.9°  40.1°  112.5°  94.9°  180°  179.6°  179.7°  156.7°  144.6°  47.2°  40.3°  180°  178.8° 

A('-13) 

165°  88.5°  138.8°  144.6°  40.3°  178.1°  179.2°  40.1°  51.5°  1 13.0°  81.4°  179.8°  42.8° 

(97) 
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ax 

"—133.1°” 

a2 

-69.3° 

a3 

99.2° 

_a4_ 

62.6° 

It  is  important  to  note  that  in  the  case  of  determination  of  all  16  elements  of  Mueller  matrix 
by  approach  Eq.(91)-  Eq.(96)  (i.e.,  k=  16),  the  absolute  minimum  of  condition  number  is 
VK  .  =  20  «  4.472 .  It  is  reached  for  condition  Eq.(76)  rewritten  in  form: 

*>16x16  >mm  ”  V  J 


Ao->«) 

a(2'-16) 

a<'-16) 


3.2°  3.2°  3.2°  3.2°  162.5°  162.5°  162.5°  162.5°  32.6°  32.6°  32.6°  32.6°  167.0°  167.0°  167.0°  167.0° 

131.1°  131.1°  131.1°  131.1°  53.4°  53.4°  53.4°  53.4°  1.3°  1.3°  1.3°  1.3°  159.7°  159.7°  159.7°  159.7° 

131.1°  162.5°  32.6°  167.0°  131.1°  162.5°  32.6°  167.0°  131.1°  162.5°  32.6°  167.0°  131.1°  162.5°  32.6°  167.0° 

3.2°  53.4°  1.3°  159.7°  3.2°  53.4°  1.3°  159.7°  3.2°  53.4°  1.3°  159.7°  3.2°  53.4°  1.3°  159.7° 
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(98) 


Taking  into  account  that  for  LC  transducers  A |. 4 e [30°,  180°]  we  get  VB^  6  min*  =  20.62  for 

conditions  Eq.(77)  and  Eq.(81)  rewritten  in  form  Eq.(98).  Thus,  13  elements  of  the  Mueller 
matrix  Eq.(93)  with  orthogonal  eigenpolarizations  can  be  determined  from  system  of  equations 
Eq.(94)  more  precisely  than  complete  matrix  in  general  case  since 

For  time-sequential  approach  Eqs.(94)  and  (96)  the  developed  software  Fig.  10  improves 
the  measuring  of  complete  Mueller  matrix  and  incomplete  Mueller  matrix  of  the  form  Eq.(93)  by 
breadboard  of  Mueller-polarimeter,  Fig.6,  with  optimal  parameters  Eq.(97)  and  Eq.(98). 

Before  testing  the  developed  polarimeter,  Fig.6,  we  performed  a  series  of  calibrating 
measurements  by  polarimeter  Fig.  5  to  determine  of  real  parameters  of  LC  transducers  (their 
Mueller  matrices)  when  driving  voltage  had  been  applied  in  accordance  with  Eq.(97)  and 
Eq.(98).  Then  elements  of  vectors  R,  S  and  matrix  Bnxn  had  been  calculated  in  accordance  with 
Eq.(92).  Here  n  is  equal  to  16  for  measurement  of  complete  Mueller  matrix  and  13  for  matrix  of 
the  structure  Eq.(90),  correspondingly. 

Finally,  by  calibrated  polarimeter  we  had  performed  measurements  of  the  Mueller  matrices 
of  linear  polarizer  (prism  Glana)  and  polymer  achromatic  quarter  wave  plate  with  orientations  0°, 
45°,  90°.  Measurements  were  carried  out  with  13  intensities  in  the  case  of  incomplete  matrix 
Eq.(93)  and  with  16  intensities  in  the  case  of  complete  Mueller  matrix.  Results  of  these 
measurements  and  error  assessments  are  shown  below  in  the  tables  1  and  2.  For  comparison 
table  3  contains  tabular  Mueller  matrices  of  91°  wave  plate  and  linear  polarizer  with  P= 0  with 
orientations  0°,  45°,  90°. 


Table  1.  Results  of  measurements  of  Mueller  matrices  of  achromatic  zero  order  quarter 
wave  plate  (A«91°)  and  linear  polarizer  (P& 0)  by  13  intensities  approach. 

quarter  wave  plate 

a=0°  (8M=  1 .4%)  a=45°(5M=  1 .4%)  a=90°(5M=  1 .4%) 
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r  1.0155  0.0008  0.0001  -0.0014^ 

0.0008  1.0139  0.0012  0.0001 

0.0001  -0.0003  -0.0177  1.0132 
v-0.0014  -0.0002  -1.0138  -0.0174, 


r  1.0157  -0.0010  -0.0005  0.0000  ^ 
-0.0010  -0.0159  0.0010  -1.0139 
-0.0005  0.0002  1.0121  -0.0008 

v  0.0000  1.0140  0.0007  -0.0186, 


r  1.0149  0.0011  0.0004  0.000  P 

0.0011  1.0141  0.0018  -0.0011 

0.0004  0.0003  -0.0180  -1.0138 
v  0.0001  -0.0015  1.0131  -0.0196, 


y=0°  (5M=1.5%) 

r  1.0155  1.0145  -0.0008  -0.0005^ 

1.0145  1.0155  0.0013  0.0028 

-0.0008  0.0002  0.0001  0.0005 

v -0.0005  -0.0006  -0.0004  0.0015  , 


linear  polarizer 

y=450(SM=1.4%) 

r  1.0157  -0.0003  1.0146  0.0003  ^ 
-0.0003  0.0015  0.0008  -0.0018 

1.0146  -0.0006  1.0127  -0.0009 
v  0.0003  0.0009  -0.0006  0.0001  , 


Y=90°  (5M=1.4%) 

r  1.0147  -1.0140  0.0003  0.0001  ^ 

-1.0140  1.0131  -0.0011  -0.0001 

0.0003  0.0008  0.0016  0.0009 

v  0.0001  -0.0006  0.0010  0.0001  , 


Table  2.  Results  of  measurements  of  Mueller  matrices  of  achromatic  zero  order  quarter 
wave  plate  (A«91°)  and  pure  linear  polarizer  (P«0)  by  16  intensities  approach. 


a=0°  (5M=1.4%) 

r  1.0154  -0.0006  0.0007  -0.0004^ 

0.0007  1.0144  0.0018  0.0004 

-0.0001  -0.0017  -0.0181  1.0129 
v -0.0001  0.0004  -1.0145  -0.0175, 

Y=0°  (5M=1.5%) 

r  1.0151  1.0145  -0.0009  0.0018  ^ 

1.0149  1.0140  -0.0010  0.0007 

-0.0016  -0.0015  0.0004  0.0021 
v -0.0010  0.0004  -0.0049  -0.0014, 


quarter  wave  plate 
a=45°(5M=  1 .4%) 


f  1.0162 

0.0008 

0.0000  -0.0005^ 

-0.0001 

-0.0188 

0.0006  -1.0133 

0.0006 

0.0008 

1.0112  0.0007 

v  0.0004 

1.0133 

0.0004  -0.0173, 

linear  polarizer 

Y=45°(5M=1.4%) 

r  1.0144  0.0002  1.0142  -0.0003^ 

-0.0001  0.0005  -0.0010  0.0002 
1.0141  -0.0001  1.0140  0.0002 
v -0.0006  0.0009  -0.0013  -0.0005, 


a=90°(5M=l  .4%) 

r  1.0149  -0.0003  0.0003  0.0007^ 

-0.0010  1.0153  -0.0007  0.0005 
-0.0007  0.0016  -0.0169  -1.0134 
v -0.0002  -0.0005  1.0124  -0.0177, 

Y=90°  (5M=1.4%) 

r  1.0154  -1.0129  0.0017  -O.OOlT 
-1.0149  1.0117  -0.0007  -0.0007 
0.0020  0.0012  0.0038  -0.0019 
v -0.0007  0.0000  -0.0014  -0.0003, 


Table  3.  Tabular  Mueller  matrices  of  91°  wave  plate  and  linear  polarizer  with  P= 0  with 
orientations  0°,  45°,  90°. 


wave  plate 


a=0° 

P  0  0  0  ^ 

0  10  0 

0  0  -0.017  1 

v0  0  -1  -0.017, 

Y=0° 

P  1  0  0^ 
110  0 
0  0  0  0 
V0  0  0  0, 


a=45° 

P  0  0  0  ^ 

0  -0.017  0  -1 

0  0  10 

v0  1  0  -0.017, 

linear  polarizer 
Y=45° 

P  0  1  0N 
0  0  0  0 
10  10 
v0  0  0  0, 


a=90° 

Poo  0  N 
0  10  0 
0  0  -0.017  -1 

v0  0  1  -0.017, 

Y=90° 

r  1  -100^ 
-110  0 
0  0  0  0 
v  0  0  0  0, 
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Error  of  measurements  described  in  the  tables  above  was  calculated  as 


SM  = 


Mtab  -Mexper 


M 


■100%, 


tab  | 


(99) 


where  Mtab  -  tabular  (exact)  Mueller  matrices  calculated  for  studied  objects  with  given 
orientations;  Mexper  -  measured  matrices.  ||  ||  -  Euclidian  norm. 

All  matrices  in  tables  2  and  3  are  nonnormalized  and  are  results  of  600  times  averaging. 
The  time  of  measurement  is  about  5  sec  for  complete  Mueller  matrix  and  about  4  sec  for 
incomplete  Mueller  matrix  Eq.(93). 

Taking  into  account  the  fact  that  the  number  of  informative  matrix  elements  can  be  less 
that  13,  see  for  example  [33,34],  the  increasing  of  the  accuracy  which  is  achieved  by  proposed 
method  can  potentially  be  higher  (several  times  or  even  order).  Indeed,  let  us  consider  the 
following  cases. 

First:  the  Mueller  matrix  describing  the  objects  characterized  by  elliptical  phase  anisotropy 
and  possibly  isotropic  depolarization.  These  can  be  the  media  of  biological  nature,  see,  for 
example,  [8]  and  references  herein.  The  matrix  has  a  form: 


"hi 

0 

0 

0 

0 

mz2 

™23 

™24 

0 

Ml  2 

Mil 

Ml4 

0 

m42 

m43 

Mm 

(100) 


The  analysis  shows  that  minimum  value  of  condition  number  for  matrix  Eq.(100)  is  13.2 
or  by  one  and  a  half  times  lesser  that  one  for  complete  Mueller  matrix.  It  is  reached  for 
conditions 


'125.5° 

175.1° 

141.9° 

153.2° 

60.4° 

179.5° 

52° 

57.4° 

51.6° 

57.4°' 

160.9° 

112.8° 

166.2° 

171.3° 

178.6° 

49.1° 

43.3° 

52.3° 

179.9° 

46.2° 

A^“10) 

40.6° 

177.4° 

179.8° 

92.1° 

46° 

111.4° 

169.7° 

72.4° 

175.3° 

63° 

a(1-1°) 

.  4  J 

50.8° 

174.3° 

51.2° 

120.9° 

150.4° 

78.9° 

65.6° 

161.1° 

90.1° 

71.9°_ 

(101) 

and: 


ax 

'  37.1°  ' 

a2 

98.7° 

a3 

90.6° 

a4_ 

_-53.!°_ 

Second  example  is 


(102) 
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mn  0  0  0 

0  m22  m23  0 

0  m32  m33  0  (103) 

0  0  0  m44 

This  matrix  describes  the  media  characterized  by  circular  phase  anisotropy  (optical 
activity)  possibly  with  isotropic  depolarization  as  well.  The  minimum  value  of  condition  number 
for  matrix  Eq.(103)  is  7.5  and,  thus,  the  error  of  measurement  of  the  matrix  Eq.(103)  decreases 
up  to  2.7  times  comparing  with  the  case  of  complete  Mueller  matrix  measurements. 

The  conditions  for  the  measurements  are  follows  : 


-Ao-6)- 

'158.7° 

176° 

40.3° 

131.2° 

o 

O 

100.8°' 

a(2‘-6) 

o 

O 

178.5° 

178° 

50.1° 

81.3° 

152.2° 

Ao-6) 

57.6° 

75.4° 

72.1° 

105.6° 

174.2° 

179.2° 

A(‘-6) 

146.9° 

119.8° 

80.8° 

40.3° 

179.7° 

91.3°  _ 

and: 


ax 

i 

00 

o 

_ 1 

a2 

91.8° 

a3 

81.4° 

a4_ 

45.1° 

(104) 


(105) 


4.  Conclusions 

Summarizing  obtained  above  results  we  can  formulate  next  conclusions. 

Anisotropy  of  homogeneous  media  with  orthogonal  eigenpolarizations  has  been 
determined.  In  particular,  when  we  use  for  modeling  such  media  an  equivalence  theorem  in  the 
form  of  Eq.(10),  eigenpolarizations  will  always  be  orthogonal  if  linear  A  and  circular 
birefringence  (p  have  the  values  Eq.(55)  and  Eq.(56),  whereas,  other  four  anisotropy  parameters 
( R ,  P,y,a)  can  take  an  arbitrary  values.  One  particular  case  needs  to  be  noted  individually:  when 
the  medium  has  simultaneously  two  types  of  amplitude  anisotropy  (circular  and  linear)  only, 
there  no  exist  any  conditions  for  such  medium  to  have  orthogonal  eigenpolarizations.  This 
conclusion  follows  directly  from  relations  Eq.(52)-  Eq.(54). 

Derived  features  of  anisotropy  of  media  with  orthogonal  eigenpolarizations  leads  to  special 
symmetry  in  respective  Jones  and  Mueller  matrices  (see  Eqs,  (14)  and  (42))  and  to  additional 
interrelations  between  matrix  elements  (see  Eqs. (21),  (23),  (25),  (29),  (30),  and  (31)).  These 
relations  can  be  useful  to  simplify  an  analysis  of  properties  of  such  media  and  can  be  utilized  for 
optimization  of  measurement  procedure.  In  present  work  we  used  an  equality  between  elements 
of  first  column  and  row  of  the  Mueller  matrix  Eq.(42)  to  reduce  the  measurement  time  with 
preservation  of  accuracy  for  polarimeter  with  four  LC  transducers.  Optimum  configurations  for 
transducers  are  presented  in  Eq.(97). 

Finally,  we  have  assembled  and  tested  breadboard  of  ID  Mueller-polarimeter  with  four  LC 
transducers  (Fig.6).  Developed  software  for  control  of  the  polarimeter  operation  allows 
measuring  the  Mueller  matrices  basing  on  16,  13  and  less  intensities  approaches  and  analyzing 
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the  measured  matrices  to  obtain  depolarization  (Cloude’s  entropy  [35,  36])  and  anisotropy  (in 
accordance  with  general  equivalence  theorem  [33])  parameters.  Thus,  “operational  flexibility”  of 
proposed  polarimeter  which  is  achieved  by  using  the  polarization  transducers  with  four  LC  and 
taking  into  account  the  exact  informativity  of  the  matrix  elements,  allows  considerably 
increasing  the  measurement  accuracy. 
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List  of  Symbols,  Abbreviations,  and  Acronyms 

(j)mn  -  phase  of  Jones  matrix  elements  tmn; 

<f>i(mn±  ki)  -  phase  of  sum(difference)  of  Jones  matrix  elements  i(tmn±tki); 

X  -  complex  variable; 

LA  -  linear  amplitude  anisotropy; 

CA  -  circular  amplitude  anisotropy; 

LP  -  linear  phase  anisotropy; 

CP  -  circular  phase  anisotropy; 

P  -  value  of  linear  amplitude  anisotropy; 
y—  azimuth  of  linear  amplitude  anisotropy; 

R  -  value  of  circular  amplitude  anisotropy; 

A  -  value  of  phase  amplitude  anisotropy; 
a-  azimuth  of  linear  phase  anisotropy; 

<p-  value  of  circular  phase  anisotropy; 

a,  b  -  real,  imagine  part  of  some  sums/differences  of  Jones  matrix  elements; 
LCVR  -  Liquid  Crystal  Variable  Retarder 
E  -  Jones  vector  of  radiation; 

S,  R  -  Stokes  vectors; 

T,  tm„  -  Jones  matrix,  Jones  matrix  elements; 

M,  mmn  -  Mueller  matrix,  Mueller  matrix  elements; 

I-  intensity  of  radiation; 

I  -  vector  of  intensities; 

So  -an  exact  Stokes  vector; 

SK  -  relative  error  of  vector  X; 

B  -  characteristic  matrix  of  stokes-polarimeter; 

G  -  characteristic  matrix  of  Mueller-polarimeter; 

V—  condition  number; 

C,Cmn  -  coherency  matrix; 

\|/  -  eigenvector  of  coherency  matrix; 

-  nondepolarizing  (deterministic)  Muller  matricx; 

jUk  -  eigenvalues  of  coherence  matrix 
□  S  -  mean  value  of  error  for  each  Stokes  parameter 
JSIV  -  peak  value  of  error  for  each  Stokes  parameter 
H-  entropy; 

9a,(p)  =  azimuth  of  analyzer,  polarizer 
P  -  azimuth  of  polarization  ellipse 
s  -  ellipticity  angle  of  polarization  ellipse 
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